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Abstract 

We consider weak solutions u £ uo + Wq' (fi, l w ) n L°° (Q,M. N ) of second order nonlinear elliptic 
systems of the type 

— div a( • , u, Du) — fe( ■ , it, Du) in Q, 

with an inhomogeneity satisfying a natural growth condition. In dimensions n £ {2, 3, 4} we show that 
M n_1 -almost every boundary point is a regular point for Du, provided that the boundary data and the 
coefficients are sufficiently smooth. 
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1 Introduction 

In this paper we are concerned with the existence of regular boundary points for the gradient of bounded, 
vector-valued weak solutions u £ W 1,2 (f2, M. N ) n L°°(Cl,WL N ) of nonlinear, inhomogeneous elliptic systems 
of the form 

— diva( • , u, Du) = b( - ,u, Du) in fi (1-1) 

with boundary values uo on dfl in the sense of traces. Here C K n is a bounded domain of class C 1:U 
and it G C 1 ' Q (n,R JV ) for some a G (0, 1). The coefficients a : £1 x R N x R nN -» R nN are assumed to 
be Holder continuous with exponent a with respect to the first two variables and of class C 1 in the last 
variable, satisfying a standard quadratic growth condition. Furthermore, we assume that the right-hand side 
b : £1 x x M. nN — »■ satisfies a natural growth condition and that an additional smallness condition on 
||u||j,oo holds. In general we cannot expect a weak solution to a nonlinear elliptic system - in contrast to weak 
solutions to a single equation - to be a classical one of class C 2 , see 115, .24] . Nevertheless, a partial regularity 
result still holds true which can be stated as follows: every weak solution u G W 1,2 (Cl, M. N ) n L°°(f2, M") 
to the inhomogeneous system is of class C 1 near a point Xq if and only if a certain excess quantity is 
sufficiently small and the mean values of u and of Du on balls B p (x ) do not not diverge for p \, 0. To be 
more precise, if we denote by 

Reg Du (Q) := {x Q e H : Du G C°(Ti n A, W lN ) for some neighborhood A of x } 

the set of regular points for Du (in the interior and at the boundary), and by Sing£) M (f2) := fi \ Heg Du (fl) 
the set of singular points of Du, then the singular set is characterized via Sing £)tl (ri) = EUE„ with 

S := \x e Q: liminf f \Du - (£»w)nns p (x ) | 2 > or limsup \{Du)nnB Jx )\ = oo\ , 

p^0+ JOnB p (xo) p^0+ J 



\x efl: limsup | {u)n nBp (x ) \ = 00 f 
1 p^o+ ' J 
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Furthermore, the gradient Du of the weak solution is locally Holder continuous with exponent a in a (small) 
neighborhood of every point xq & Reg Du (Q), see [55] (and [30] IS] f° r the non-quadratic analogues). This is 
the up-to-the-boundary extension of the interior partial regularity results obtained in various papers starting 
from [23] [21] [31]. By Lebesgue's differentiation theorem, the regularity criterion stated above applies to 
almost every point in fl, whence |ri\Reg £3u (r2)| = 0. However, this does not yield the existence of even one 
single regular boundary point for weak solutions to general nonlinear elliptic systems, since the boundary 
dQ itself is a set of Lebesgue measure zero. By contrast, due to Giaquinta's counterexample |20j . it is well 
known that singularities may occur at the boundary even if the boundary data is smooth. 

The question of dimension reduction of the singular set (in the sense that it is not only negligible with 
respect to the Lebesgue measure but that its Hausdorff dimension is bounded strictly below n) has received 
considerable attention in recent years. Some significant results were first obtained for weak solutions to 
systems satisfying special structure conditions: for quasilinear systems of the form 



various partial regularity results were established, stating that the weak solution u (instead of its first 
derivative) is locally Holder continuous. To bound the Hausdorff dimension of the singular set Sing u (f2), 
we recall that a regular point xq G of u is a point where u is locally continuous and is characterized via 
a smallness condition on the lower order excess functional 



e. g., see [23j [TH [38] [25] [3J . Since the set of non-Lebesgue points of every H /1,p -map has Hausdorff dimension 
not larger than n — p, the Hausdorff dimension of Sing u (£l) cannot exceed n — 2. If the coefficient matrix 
a(-, •) of the quasilinear system is further assumed to be of diagonal form, it is known that the weak solution 
is a classical solution (see [42] where boundary regularity is included). Useful estimates for the singular 
set are also available for nonlinear elliptic systems obeying special structure assumptions: for instance, 
Uhlenbeck established in her fundamental paper [41] a strong maximum principle for the gradient Du of 
weak solutions to nonlinear systems, provided that the nonlinear part of the coefficient function only depends 
on the modulus of Du. This was the key to obtain everywhere-regularity for Du. For an extension to the 
nonquadratic case we refer to |40[ [T] . However, neither could Uhlenbeck's techniques be carried over to the 
boundary, nor is a suitable counterexample available in the literature, leaving the question of full boundary 
regularity open for such systems. Turning the attention to general nonlinear elliptic systems, we observe 
that a direct comparison technique allows to infer local Holder continuity of the weak solution outside a set 
of Hausdorff dimension n — p in low dimensions n < p + 2, see [lOi 1111 [4] [5] [8]. By contrast, in arbitrary 
dimensions n the reduction of the Hausdorff dimension of the singular set Sing £)u (f2) for the gradient Du 
was a long-standing problem. It was finally tackled by Mingione [37]: he studied the interior singular set 
Sing£, u (f2) in the superquadratic case p > 2 for systems without ^-dependencies and with inhomogeneities 
obeying a controllable growth condition, and he succeeded in showing that the Hausdorff dimension of 
Sing £)u (51) is not larger than n — 2a. In [36] he extended these results to systems with inhomogeneities 
under a natural growth condition, covering also systems explicitly depending on u, provided that n < p + 2 
is satisfied. 

We now return to the existence of regular boundary points: we first observe that for this aim the almost- 
everywhere regularity result has to be improved to a bound for the Hausdorff dimension less than n — 1 
because this yields immediately that almost every boundary point is regular. Consequently, our objective 
is to identify additional assumptions on the coefficients or on the space dimension which guarantee this 
dimension reduction. A result in this direction was recently obtained by Duzaar, Kristensen and Mingione 
[18] : they considered weak solutions u £ W 1,p (ft, M. N ), p e (l,oo), of the homogeneous Dirichlet problem 
corresponding to (|l.ip and developed a technique which allows to carry the estimates in [37J up to the 
boundary, implying in particular the existence of regular boundary points, provided that n — 2a < n — 1 
(or equivalently a > \) is satisfied. More precisely, the authors obtained for every a 6 (5,1] that almost 
every boundary point is regular if the coefficients a(x, z) have no it-dependency or if n < p + 2 holds. In the 
quadratic case they improved this result in two ways: on the one hand, inhomogeneities with controllable 



div (a( • , u) Du) = b( - ,u, Du) , 
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growth were included, and on the other hand the condition on a was sharpened to a > i — e for some 
number e > stemming from an application of Gehring's lemma. We further mention that various results 
establishing better estimates for the (interior) singular set of minimizers of variational integral can be found 
in [33j EH- 

The main result in this paper is an extension of the result [18j to bounded weak solutions to inhomoge- 
neous systems under a critical growth condition on the inhomogeneity (giving also an alternative proof of 
[T51 Theorem 1.3]), namely the improvement of the estimate | Sing £)tl (ri)| = in the following sense: 

Theorem 1.1: Consider n £ {2,3,4} and a > i. Let f2 C R" be a domain of class C 1,a and u £ 
C 1 ^^!*). Assume that u G u + Wq' 2 (CI,R n ) n L°°(Cl,R N ) is a weak solution of the Dirichlet problem 
under the assumptions (Hl)-(HS) and (B) from Section^ and suppose that HitHioo/n.R-") < M for 
some M > such that 2L 2 M < v. Then IK™" 1 -almost every boundary point is a regular point for Du. 

This result was presented as a part of the author's PhD thesis [7J where most of the proofs and calculations 
are discussed in detail. In addition, some extensions and open questions concerning the dimension reduction 
of the singular set are collected in Section [51 In particular, in view of an observation by Kristensen and 
Mingione [33] , it is possible to replace a part of the Holder continuity assumption on the coefficients with 
respect to the x-variable (in the sense that it is only required with an arbitrary exponent) by an additional 
fractional differentiability assumption on the map x <-> a(x,u,z), see Theorem 16.41 

We close this introductory part with some remarks about the ideas behind the arguments and the 
techniques used within this paper. The strategy can be described as follows: To simplify matters we 
initially consider coefficients of the form a(x,z): If they are Holder continuous in x with arbitrarily small 
exponent, we know dimjf (Sing £) , tl (f2)) < n. If they are instead Lipschitz-continuous, then standard difference 
quotients reveal Du £ W 1 ' 2 ^, M. nN ) which implies that dimjc(Sing £)u (n)) < n — 2. Therefore, the upper 
bound on the Hausdorff dimension of Sing Dn (f2) reflects the regularity of the coefficients in x. This gives the 
impression that the regularity of the coefficients is related not only to the regularity of the solution (namely 
the local Holder continuity of Du to the same exponent), but also to the size of the singular set. Working 
from this observation, Mingione 37, 36 introduced a remarkable new technique and accomplished in the 
interior an interpolation between Lipschitz continuity on the one hand and Holder continuity on the other: 
for general a-H61der continuous coefficients the existence of higher order derivatives of the weak solution 
cannot be expected, but it is still possible to differentiate the system (|1.1|) in a fractional sense. This leads 
to the desired upper bound n — 2a for the Hausdorff dimension. If the coefficients a(x, u, z) now depend 
explicitly on u, the situation becomes more complex and the estimates are technically much more involved. 
To follow the line of arguments above we have to investigate the regularity of the map x i— > (x,u(x)). If 
the weak solution it is a priori known to be everywhere Holder continuous then x i-> (x, u(x)) is also Holder 
continuous and the arguments apply with only marginal modifications. However, in general this map is no 
longer continuous, because u may exhibit irregularities. Nevertheless, at least in low dimensions n < p + 2, 
local Holder continuity of weak solutions is guaranteed outside of closed subsets of Hausdorff dimension less 
than n — p. In other words, the set of points where u is not continuous - and where x i— ¥ (x,it(x)) is not 
regular - has sufficiently small Hausdorff dimension, hence, restricting the analysis of Du to the regular set 
Reg u (f2) of u, we still arrive at a good result for dim^Sing^ffi)), see Theorem 16. II 

In the interior this method relies essentially upon finite difference operators, fractional differentiability 
estimates for the gradient Du and interpolation techniques dating back to Campanato |12[ [9] , combined in a 
delicate iteration scheme (applied for elliptic and parabolic systems [371 (Ml HE] ) , and the necessary estimates 
are deduced by testing with (differences of) the solution. At the boundary some severe problems are caused 
by the fact that testing is allowed only for differences in tangential direction: hence, the normal direction 
still has to be recovered by exploiting the system of equations (which follows immediately if second-order 
derivatives exist). This problem was overcome first for homogeneous elliptic system (and inhomogeneous 
systems under controllable growth) by an indirect approach developed by Duzaar, Kristensen and Mingione 
[18) : via a regularization procedure involving both the original coefficients a(-, •, •) and the specific solution 
u, a family of comparison maps is constructed for which the existence of second-order derivatives is known. 
This allows to gain higher integrability for Du which in turn is used to improve the integrability of the 
comparison map by means of Calderon-Zygmund estimates (provided in [33J [Z] ) in the next iterative step. 
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When trying to apply this approach for inhomogeneous systems under critical growth, several critical 
difficulties arise: most importantly the propagation of higher integrability via the Calderon-Zygmund theory 
seems not to be clear since the natural growth condition merely gives L 1+s for the right-hand side with some 
(small) 8 > (coming from the higher integrability of Du) rather than the necessary prerequisite L?/(p _1 ) for 
some q > p. For this reason we exploit the system differently and replace the indirect comparison principle 
by a direct method, introduced by Kronz [35j as a promising approach for up-to-the-boundary regularity 
results including upper bounds for the Hausdorff dimension of the singular set, with the flexibility to attack 
higher order systems. Kronz observed that estimates for the tangential differences suffice to control the 
averaged mean deviation with respect to mean values taken over slices in tangential direction. Using an 
alternative definition of fractional Sobolev spaces based on pointwise inequalities, this helps to deduce a 
fractional differentiability property for the system, which then gives further information on the gradient 
of the solution. The overall strategy remains unchanged, i. e. existence of regular boundary points is 
still proved by a dimension reduction argument for the singular set Sing^^O): The key tool here is the 
observation that if Du belongs to a fractional Sobolev space W ' p , then the characterization of Sing £)u (0) 
and a measure density result allow to conclude that the Hausdorff dimension of Smg£>„(0) does not exceed 
n — Op. The proof of such a fractional differentiability estimate for Du is now sketched in a series of steps: 

Strategy of the proof: 

Simplifications: It suffices to consider the model situation f2 = B + and solutions u 6 W 12 (B + , R N ) n 
L°°(B + ,M. N ) which vanish on the flat part of the boundary. The general situation then follows from a 
transformation argument. Furthermore, we assume Holder continuity of u on B + with some exponent 
A > 0. This is justified by the fact that the solution is Holder continuous outside a set of Hausdorff 
dimension n — 2 in dimensions n € {2, 3, 4}. 

Tangential differences: Testing the system with differences of the solution up to the boundary is only 
allowed for tangential directions (because zero boundary values on the flat part are maintained for the test 
function). Taking into account the assumptions on the coefficients and the inhomogeneity, we end up with 
an integral estimate for \Du(x + he s ) — Du(x)\ for all unit directions e s _L e n , telling that its L 2 -norm 
decays like c\h\ aX / 2 (with a denoting the Holder exponent of the continuity condition on the coefficients 
with respect to the first and the second variable). 

An estimate for tangential derivatives: If finite differences of the full derivative Du are estimated, then 
it is reasonable that also normal differences of only the tangential derivative denoted by D'u are estimated 
similarly (if we think of Lipschitz-continuous coefficients a(x, Du) for example, this observation is trivial 
since the previous step yields the existence of second order derivatives D' Du = DD'u). This is in fact true 
(up to a small loss in the power of \h\) 7 and we thus get a first fractional differentiability estimate for D'u. 

Towards the normal derivative: Information about D n u can only be gained out of the system (in case 
of Lipschitz-continuous coefficients a(x,Du), the existence of the second order normal derivative D n Du is 
obtained from the system of equations —D n a n (x, Du) = J2k=i Dkak(x, Du)+b(x, u, Du) in a standard way). 
Looking at the simple example — div (f(x,u)Du) = b(x,u, Du) we get a first idea on how the coefficients 
might serve to improve the differentiability of D n u, because we then have a n (x,u, Du) — f(x,u)D n u, 
meaning that a n (x,u, Du) is the missing normal derivative up to a Holder continuous perturbation (a 
similar property holds true for the general coefficients). For the moment let us concentrate on a n (x, u, Du): 
mimicking the differentiable situation to a certain extent, we show by means of the estimates for tangential 
differences of Du that slice- wise mean values of a n (x,u(x), Du(x)) are differentiable in the weak sense in 
the ^-direction, and as a consequence, we obtain that the map x i-> a n (x,u(x), Du(x)) is in a fractional 
Sobolev space. 

An estimate for the normal derivative: Taking advantage of the ellipticity and the boundedness condition 
assumed for the coefficients, we find that differences of D n u are essentially dominated by those of a n (x, u, Du) 
and of the tangential derivative D'u. Together with a corresponding estimate for the tangential derivatives 
of u, this leads to a fractional differentiability result for the full gradient Du e W " 1 *' 2 for every 7 < 1. 

Getting rid of X: By an interpolation technique we gain higher integrability out of the fractional differ- 
entiability of Du. This in turn is used to improve the differentiability of Du in a suitable iteration procedure 
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up to the final result Du £ W a1 - 2 for every 7 < 1. 
2 Structure conditions and notation 

We impose on the coefficients a : f2 x x W lN — > M. nN standard conditions (here stated for general p- 
growth, even if we will concentrate on case p = 2): the mapping z M> a(x, u, z) is a continuous vector field, 
and for fixed numbers < v < L, p £ (1, 00) and all x, x £ 0, u, u £ R N , z £ R nN , the following growth, 
ellipticity and continuity assumptions hold: 

(HI) a has polynomial growth and is differentiable in z with continuous, bounded derivatives: 

\a(x,u,z)\ + (l + \z\ 2 )^ \D z a(x,u,z)\ < i (l + ^l 2 )^ , 
(H2) a is uniformly strongly elliptic, i.e. 

D z a(x, u,z) X - X > v(l + \z\ 2 ) E ^\X\ 2 VAeR nAr , 

(H3) There exists a nondecreasing, concave modulus of continuity ui a : R + — > [0, 1] 
such that co a (s) < min{l, s a } for all s £ R + and 

p-i 

\a(x, u, z) — a(x, u, z)\ < L (l + |z| 2 ) 2 cj a (\x — x\ + \u — u|) . 

The latter condition (H3) prescribes uniform Holder continuity with respect to the (x, w)-variable with 
Holder exponent a (for fixed z). Moreover, we assume the inhomogeneity b: fl x ~R N x R nN — > to be a 
Caratheodory map, that is, it is continuous with respect to (u, z) and measurable with respect to x, and to 
satisfy a natural growth condition of the form 

(B) there exists a constant L 2 (possibly depending on M > 0) such that 
\b(x,u,z)\ <L + L 2 \z\p 
for all x £ n, u £ R N with |u| < M, and z £ R nN . 

We further make some remarks on the notation used below: 

(Half-)Balls, cubes and cylinders: We write B p (y) — {x £ R™ : \x - y\ < p] and B+(y) = {x £ W 1 : 
x n > 0, \x — y\ < p} for an n-dimcnsional ball or the intersection of the ball with the upper half-space 
M™ -1 x E + , centered at a point y £ M™ (respectively £ R™ -1 x Rq in the latter case) with radius p > 0. In 
the case y = we set B p := B p (0), B := B\ as well as B+ := B+(0), B + := B^ . Furthermore, we denote 
by D p (y') the (n - l)-dimcnsional ball D p (y') := {x £ R"" 1 : \y' - x'\ < p} for y' £ R"" 1 , and by Z p (y) 
the open cylinder on the upper half-plane R™ _1 x R+ 

Z p {y) ■= D p {y') x (max{0,y„-p},j/„ + p) =: D p (y') x I p (y n ) 

for a center y =: (y 1 , y n ) £ W 1 with y n > 0. Similarly as for balls, cubes with center y £ R™ and side-length 
2p are denoted by Q P (y), upper half-cubes by Q+(y), and we further write Q° p (y) = dQ+(y) n R"" 1 x {0} 
(with the corresponding abbreviations for y = and if p = 1). 

Function spaces: We will work with functions belonging to the Holder spaces C 1 '", a £ (0, 1), and the 
(fractional) Sobolev space W e - p , 9 £ (0, l],p £ [l,oo). The definition for noninteger values of 9 and some 
preliminary material is collected in the next section. Moreover, we introduce the following notation for 
W 1,p -functions defined on a upper half-cubes Qt(y) which vanish on Q ( p (y) (in the sense of traces): 

W^(Qj(y),R N ) := {u £ W 1 * (Q+(y), R N ) : u = on Q° p (y)} . 

where y n < p is satisfied; the subspace of functions vanishing on the whole boundary is denoted by Wj' p . 
Sometimes, it will be convenient to treat the tangential derivative D'u := {D\u, . . . , D n _\u) and the normal 
derivative D n u of a Sobolev function u separately. 
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Measures and mean values: For a given set X C K fe we write L k (X) = \X\ and dimj{(X) for its k- 
dimensional Lebesgue-measure and its Hausdorff dimension, respectively. Furthermore, if h £ L X (X, M. N ) 
and < \X\ < oo, we denote the average of h by (h)x = j- x hdx, and when working on cylinders we will 
use the abbreviation (v) Xo _ p := {v)z p ( Xo )- We further define the slice-wise mean value of u in D r ((xo)') at 
almost every height x n £ I p ((xo) n ) via 

(v) x > p (x n ) := f v(x',x n )dx' . 

JD p {{x )>) 

The constants c appearing in the different estimates will all be chosen greater than or equal to 1, and 
they may vary from line to line. 

3 Fractional Sobolev spaces and finite differences 

In what follows, we will use the notation of [2] (see also [33J [18]). For a bounded open set A C K™, 
parameters 9 £ (0,1) and q £ [1, oo) we write u € W 9 ' q (A, R N ) provided that u £ L q (A,R N ) and the 
following Gagliardo-type norm of u dehned as 

Hx)\ q dxy + { J a J a \ x _ y \ n+q9 dxdyj 

is finite. In order to formulate a general criterion for a function to belong to a fractional Sobolev space we 
introduce the finite difference operator r e) h with respect to a direction e G B\ C M" and with stepsize h £ M. 
via 

T e , h G(x) ee r e , h (G)(x) := G(a: + /ie) - G(x) 

for a vector valued function G : A — > H N (this makes sense whenever x, x+he £ A). If e = e s , s G {1, . . . , n}, 
is a standard basis vector, we use the abbreviation r S: h instead of T es: h- These finite differences are related 
to the fractional Sobolev spaces (in the interior as well as in an up-to-the-boundary version) via the next 
lemma: 

Lemma 3.1 ([32]. Lemma 2.5; [15] . Lemma 2.2): Let G £ L q (Q^,R N ), q > 1, and assume that for 
9 G (0, 1], M > and some < r < R we have 

n 

V/ \r Sih G\ q dx < M q \h\ qe 
s =i J Qt 

for every h £ M. satisfying < \h\ < d where < d < min{l,i? — r} is a fixed number. In the case s = n 
we only allow positive values of h. Then G £ W b ' q (Q^ , R N ) for every b £ (0,9) and p < r. Moreover, there 
exists a constant c = c(n,q) (in particular, independent of M and G) such that the following inequality 
holds true: 

[ f ^- G ^ q dxd y< C ( Mq£qi " b) +ML f \ G \ q dx ] 

where e := min{r — p, d}. In the interior the same result holds true without any constraint on the sign of 
h with respect to the direction of the differences t S) /j. Moreover, we can consider (half-)balls or cylinders 
instead of cubes. 

In the case where G is the weak derivative of a W 1,q function v and where an estimate for finite differences 
only in tangential direction is known, we are still in a position to state a fractional differentiability result 
which is limited to the tangential derivative of v: 

Lemma 3.2: Let v £ W 1 ' q {Q'j i , R N ), q > 1, and assume that for 9 £ (0, 1], M > and some < r < R 
we have 
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for every h G M satisfying < \h\ < d where < d < min{l,i? — r} is a fixed number. Then D'v — 
(Div, . . . , G ^''(Q+.M^" 1 ^) /or every b G {0,9) and p<r. 

Proof: We first fix 6 G (0, 9) and p G (0, r). We consider arbitrary numbers h! G M + and feel satisfying 
< \h\, \h'\ < min{<i, ^-jp}. Then, using Young's inequality, standard properties of the difference operator 
and the assumption (|3.1[) on finite differences in tangential direction, we conclude for every e G (0, 9) and 
s G {l,...,n-l}: 

/" \ TnhlTshTs _ hV \1 d X 



\r s .hT s .-hv\ q dx 

+ 



< (\ti\-*\h\- e * + \h\-*- e *) / \T n , h ,r s , h r s< - h v\ q dx 

< 2 \h'\~ 9 \h\- eq [ \r s . h T n . h ,v\" dx + 2 \h\- q - 6q [ 

J Qt- d J Q,- 

< 2\h\- 0q f \T sh D n v\ q dx + 2\h\- 0q [ \r sh D s v\ q dx < 4M' 
uniformly in h, h! . From |16l Lemma 2.2.1] we infer (for possibly smaller values of \h 

\ h l\-{e-e)q 1^1-9 f \ Tnh ,T sh v\ q dx < c( [ \Dv\ q dx + M q 
Jot K JQt 



r-2d 



and the constant c depends only on 9, q, e, d and r — p. Considering the limit h — > 0, we hence end up with 
\h'\-(o-e) q f \ Tnh ,D s v\ q dx < c( [ \Dv\ q dx + M q 

Keeping in mind that the index s G {1, . . . , n — 1} is arbitrary, we combine the latter inequality with (|3.1 
to find 

V/ \r sh D'v\ q dx < c\h\ {e - £ ^ q ( [ \Dv\ q dx + M q 



S=l " V T ._ 2 d 

1 p 



for all h G M satisfying < \h\ < min{d, ^-3-^} where we only allow positive values of h if s = n. For 
e = (9 — b)/2 the application of Lemma I37T1 with 9, r replaced by 9 — e, r — 2c? finishes the proof. □ 

The following interpolation inequality can be found in [9l Lemma 2.V] and is essentially based on the 
inequality in [12l Theorem 2.1] for the case p = 2. 

Theorem 3.3: Let X,9 G (0,1], p G (l,oo) and u G C°' X (Q,R N ) such that Du G W e > p (Q, R nN ) with 
p9 < n, where Q C is an (upper) cube. Then 

Du£L s (Q,W lN ) for all s < + „? . 

n — p9\ 

Moreover, 

J \Du\ s dx < c(n, N,p,9, A, s, \Q\, |MIwi+9.i>(q,rW), Mc°.*(q,r^)) • 

The next lemma enables us to conclude from difference estimates for a map v an appropriate estimate 
for the averaged mean deviation with respect to slice-wise mean values: 

Lemma 3.4 ([35]): Let a < -|, n > 2, r > 0, Z p (xq) C Q + for some xq G Q + UQ°. Furthermore, assume 
that v G L p (Z p (x ),R N ), p > 1, satisfies 



I \r h .ev\ p dx < KP\h\ T P 

J Z ap (x ) 
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for some K > 0, all e G S n 1 with e _L e n and ftgE ?ot£/i \h\ < 2ap. Then, for every /3 G (0, t) there exists 
a function F G L p (Z a p(x$)) such that 

\F\ p dx < c{n,p,T,j3)K p p^-fiP 

' Z ap (x a ) 

and 

l I Hx',x n )-v(y',x n )fdy'dxY<c(n,f3)r^F(z) 

>Z r (z) JD r (z>) ' 

for every exponent p G [l,p), almost all z G Q + U Q° and all r > such that Z r [z) C Z ap (xo). 

A different definition for fractionai Soboiev spaces, based on pointwise inequaiities, can be derived as 
follows: Let O C R™ be a bounded domain, p > 1 and £ (0, 1]. Following the approach of Hajlasz in [28], 
we set 

/) := {. 9 G i p (fi) : 3E 1 C f2, = such that 

- /(y)| < |x - + g{y)) for all a;, y e il \ E) , 

and we define the fractional Soboiev space via 

M e ' p (Q,M N ) := {/ e L P (Q,R N ): D e ' v {VL-f) ^ 0} . 

M 8 >P(n,R N ) is equipped with the norm 

ll/llMe.p(o,R«) : = ll/llLp(n,R«) + pf s Hs'lUpfn) ■ 

g6D e .p(f2;/) 

For p G (l,oo), due to the convexity of L p , to every / G M e ' p (Cl, WL N ) there exists a unique function 
g G L p (£l) which minimizes the L p (f2)-norm amongst all functions in D ' p (f2;/). We highlight that this 
definition has its origin in the definition of Soboiev spaces in the context of arbitrary metric spaces (replacing 
\x — y\ by dist(x, y)) and that it does not use of the notion of derivatives (for a more detailed discussion 
of the metric setting we refer to |29j). Employing the Hardy-Littlewood maximal function we see that 
this "metric" Soboiev space coincides with the classical Soboiev space for the integer order 9 = 1 and 
sufficiently regular domains (e.g. with Lipschitz boundary). More precisely, provided that p > 1, there 
holds M 1,p (fi,R ) = W 1,p (tt, R. N ) for all bounded domains ft with the so-called extension property, 
meaning that there exists a bounded linear operator E : W 1,p (Ci, R ) — > W 1,P (W\ R ) such that for every 
/ G VF 1,p (r2, R N ) there holds Ef = f almost everywhere in ft. Instead, the equivalence fails if p = 1, see 
[2 7) . Furthermore, the definitions of the classical and the metric fractional Soboiev spaces immediately yield 
for all bounded domains ft, fractional orders 9 G (0, 1) and p G [1, oo) the following inclusion: 

M e ' p {il, R N ) C W e '' p {n, R N ) for all 9' G (0, 9). 

The following lemma provides an integral characterization of fractional Soboiev spaces for domains 
satisfying the mild Ahlfors regularity condition, which demands the existence of a positive constant fco such 
that 

(Kq) \B p (xq) n il\ > fcn p n f° r a U points xq G £1 and every radius p < diam(J7) . 

In other words: the domain is not allowed to have external cusps. We note that the latter condition is for 
example satisfied by the large class of domains with Lipschitz-continuous boundary. 

Lemma 3.5: Let Q C R™ be a domain which fulfills an Ahlfors condition (Kq), 9 G (0,1], p G (l,oo). 
Then the following two statements are equivalent: 

(i) f g M°'P(Q,R N ) 
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(ii) f G L 1 (f2,R Ar ) and there exists a function h G L p {Vl) and a radius Rq > such that 

/ l/-(/K(*„)nn|<k < AfoO (3.2) 

JB p (x a )nn 

for almost all xq G and p < i?o • 

Proof: The implication (i) =>• (ii) follows by standard properties of the Hardy-Littlewood maximal function 
for the choice h = AM{g) with g G £) e,p (f2; /). The reverse implication (ii) =>• (i) is an easy adaptation of 
the proof of Campanato's integral characterization of Holder continuous functions, see e. g. [5H1 Chapt. 1.1, 
Lemma 1]. □ 

Remarks 3.6: In fact, the following local version of the integral characterization holds: let Xq G fl and 
R > such that 



L 



1/ - (/)s r («)nn| dx < r K z ) 

'B r (z)nn 

for almost all z G fi, B r (z) C Bi?(a;o) and h G U>{il) as above. Then there holds / e Af e ^(B fl . /2 (a;o)n^, l w ) 
with 

~ f(y)\ < c(n, ha, 0) \x~y\ e (h(x) + h{y)) 

for almost all x, y G Bri^xq) n fi. In view of Jensen's inequality and the fact that the Hardy Littlewood 
maximal operator is a bounded map from L p to itself, this characterization allows to infer the inclusion 

W e ' p (n,R N ) C M e,p (Q,M JV ) 

whenever J7 satisfies an Ahlfors condition (Kn), 9 G (0, 1) and pG (1, oo). 

Moreover, we note that (i) implies indeed the following statement: there exists a function h £ L p (f2) 
and a radius i?o > such that 



for all q < p and almost all Xq £ f2 and p < Rq 



(/ l/-(/)B,(« )nnl 9 d*)* < /M^o) 



4 Some basic facts about the solution 



In what follows, we restrict ourselves to the model case fl = Q2 , and we study weak solutions u G 
W^ p (Qt,R N )r)L QC (Q+,R N ) of the system 

— diva( ■ , u, Du) = b(-,u,Du) in Q\ . (4.1) 

By a transformation argument this covers the situation of general inhomogeneous systems of type (|1 . 1|) on 
arbitrary domains ft of class C 1,a . Moreover, we argue under the permanent assumption that the weak 
solution u of system (|4.1[) is Holder continuous on Q + with Holder exponent A for some A G (0, 1). This 
assumption will later be justified by the fact that in low dimensions the weak solution u is a priori known to 
be Holder continuous outside a set of Hausdorff dimension n — 2 (and since we are interested in the behavior 
of Du on the boundary which is of Hausdorff dimension n — 1 this information is sufficient to forget about 
the bad set where u is not Holder continuous). 

We now present some tools needed in the remainder of the paper: first, we recall the well-known 
Caccioppoli inequality in an up-to-the-boundary version. The fact that the oscillations of u are due to its 
continuity arbitrarily small in a cylinder - provided that the side length of the cylinder is chosen sufficiently 
small - allows to simplify the estimates which are usually slightly more involved for nonlinear elliptic 
systems with inhomogencities under a natural growth condition. As a matter of fact we here do not need 
the smallness assumption |u| < M with 2L 2 M < v. 
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Lemma 4.1 (Caccioppoli inequality revised): Let u G W r (Qjj,R ) n£ 0o (Q2>"* ) ^ e a wea ^ solu- 
tion of (|4.ip under the assumptions (H1)-(H3) and (B). Assume further u G C°' A (Q + , 1*). XTien £/iere 
errasi positive constants c cacc = c cacc (n, N,p, ^) and p cacc = p caC c(p, ~, -7?: ^> Mc°> a (q+,r jv )) suc/l that 
for every £ G and every cylinder Z p (y) C Q + im£/i y G Q + U Q° and y„ < p < p ca cc there holds: 

/ - 7(C ® e„)| 2 dx < ? CQCC ( / 2 dx + p 2Q (1 + |e|) P+2a ) . 

Jz p/2 (y) K JZ p (y) V P ' ' 

Here we have used the T^-function which is in general defined by V(£) = (1 + \{,\ 2 Y P ~ 2 ^ 4 £, for all £ G R fc 
for some k G N (in the quadratic case it is just the identity map) and which is in particular a bi-Lipschitz 
bijection on R fc . Secondly, we recall an estimate concerning finite tangential differences of Du which is the 
starting point to proceed to fractional differentiability estimates for Du and hence to dimension reduction 
arguments for the singular set: We consider 8 G (0, 1) and assume u G W^ ,P (Q2, R N ) n L°°((3j, WL N ) to be 
a weak solution of system (|4.1j) . Then for every cut-off function 77 G C{j°(Qi-s, [0, 1]) and every tangential 
direction e G S 1 " -1 with e _!_ e n there holds 



/ n 2 \r e . h V(Du)\ 2 dx < c(\h\ 2a [ 



(1 + \Du(x)\ p + \Du(x + he)\ p + \h\- p \T^ h u{x)\ p ) dx 

Q+nspt(r;) 

(1 + \Du{x)\ 2 + \Du(x + he s )\ 2 )^ \r eM u{x)\ 2a dx 

-nspt(r;) 



+ j Q+ ( 1+ \ DU ( X )\ P ) \Te,-h(V 2 TeMx))\d^ 



(4.2) 



for all h G M with \h\ < 6, and the constant c depends only on n, N,p, —, ||u||loo and ||Z)^||x,oo. 
We highlight that this estimate is the up to the boundary analogue of estimate (4.7)], and its proof 
follows the line of arguments in [36] : Testing the weak formulation of (14. ip with the function T e ^hip for 
f G W / q 1:P (Q + , 1*) n L°°(Q + ,Hl N ) with sptyj C Qi-a, we first use partial integration for finite differences 
on the left-hand side which results in integrals involving r Sj / l (a(x, u(x), Du(x))) . Decomposing 

r ei h(a(x, u(x), Du{x))) 

= a(x + he, u{x + he), Du(x + he)) — a(x, u{x + he), Du{x + he)) 

+ a(x, u(x + he), Du(x + he)) — a(x, u(x), Du{x + he)) 

+ a(x, it(x), Du(x + fte)) — a(x, w(x), Du(x)) 
=: A(h)+'B(h) + e(h), (4.3) 



[,A(/i) + £(/i) + e(hj] ■ Dipdx = / b(x,u,Du)-T e - h <pdx. (4.4) 



we hence find 



Choosing ip = i] 2 T e _hU, we have to estimate the various terms by taking advantage of the growth and 
continuity assumptions of the coefficients and the inhomogeneity exactly as in [36) . and we then end up 
with the desired inequality (|4.2I) . 



5 The proof of Theorem 11.11 

5.1 Higher integrability of finite differences 

We first state a higher integrability estimate for both Du and for finite differences of Du (again motivated 
from |36)V which will allow later to end up with a slightly sharper estimate on the Hausdorff dimension 
of the singular set. We first observe the well-known existence of a higher integrability exponent sq > 2 
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depending only on n,N, ^ and [u]c°. a (q+,r jv ) such that u 6 T4^ 1,s ° (Q+ , R^) for all p < I. Furthermore, 
for every center xo G Q + U Q° and every radius p £ (0,1 — \xq\) there holds 

— /* - 

c) s ° < c(n,N,^^,[u] C o, HQ+ ^)) (f (1 + \Du\ 2 ) dx)* , (5.1) 



|£>u| so dxJ 

see e.g. [SI Lemma 4.1]. Combining the higher integrability with (|4.2I) we obtain similarly to 36, Section 5, 
step 2] a higher integrability result for T e ^.Du: 

Proposition 5.1: Let u G W^ 2 (Q^ ,R N )r\L°°(Q+ ,M. N )nC - x (Q + ,R N ) be a weak solution of flUi} under 
the assumptions (H1)-(H3) and (B). Furthermore, let Z p (xq) C Q + for some x$ G Q + U Q° , a G (0, j^), 
e G S 1 ™ -1 with e _L e„ and h G K wi£/i G (0,2crp). Then there exists a higher integrability exponent 
s G (2, So) depending only on n,N, — , and [ - u]c°. a (q+,r jv ) such that 



\T eJl Du\ s dx < c\hf¥- (/ (l + \Du\ 2 ) dx) 2 
for a constant c= c(n,N, ^f, [u] c o,x {Q + tKN) , p,a) . 

Proof: We consider in the sequel the tangential directions e G S n_1 , i.e. e _!_ e n , and we initially look 
at numbers h G K satisfying < 1. Recalling the abbreviations for A(h), 23(/i) and C(/i) from (|4.3[) . 
representing the differences of the coefficients a(-, •, •) with respect to each variable, we set 

T e , h U ~a(h\ Z^W m/n ~3(/Q 

and we define G{h) = J D z a(x,u(x),Du(x) + tT e ,hDu(xyj\ dt. Dividing the previous identity (14.4)) by 
\h\ aX ' 2 (which is half the power of | h | to be expected in (I4.2|) for A-H61der continuous solutions) we get 

6(h) Dv h - Dip dx = / [A{h) + %(%)]■ Dip dx + I \h\~^ b(x,u,Du) ■ r e - h ipdx (5.2) 

Q+ JQ+ JQ+ 

for all functions (p G W^ 2 {Q+_ W ,M. N ) D L°°(Q+_ lh{ ,R N ), i.e. the map v h G W X ' 2 (Q+_ W ,M N ) is a weak 
solution to the linear system (15.21) for every h G R with < 1. In the next step we infer Caccioppoli-type 
inequalities for the functions Vh, for which the constants may be chosen independently of the parameter h. 
For this purpose we first observe some simple properties due to (H1)-(H3) and the Holder continuity of u 
with exponent A: 

\A(h)\ < L (l + \Du{x + he)\) , 

\$(h)\ < L [v]% , HB+>MN) (1 + \Du(x + he)\) , 

^|A| 2 < e(h)\® A < L|A| 2 VAel' lJv . 

For a, p and Xq fixed according to the assumptions of the proposition, we next choose h G M such that 
\h\ G (0,2ap) and consider intersections of balls B^(y) with the upper half-plane R™ -1 x R + for centers 
y G Z( 1 _ r7)p/2 (xo) satisfying B^(y) C (implying that < R < 1 - \h\ - max fee{1) ... )n} \y k \) and 

Vn < i- e - we fi rs t study the situation for centers close to the boundary. Furthermore, we take a cut-off 
function r/ G C^ 3 (B 3R / i (y), [0, 1]) satisfying n = 1 on B R / 2 (y) and |Dry| < -|, and we choose ip := 77 2 Wft as a 
test function in (|5.2[) . Taking into account Dp = n 2 Dvh + 2rjVh <8 Dr\, we estimate the various terms arising 
in (|5.2p : using Young's inequality with e G (0, 1) and the estimates for A(h), H(h) and C(h) given above we 
see 



• v\ n 2 \Dv h \ 2 dx< / n 2 e{h)Dv h - Dv h dx, 

JB+(y) JB+(y) 

• ( 2ri\e(h)Dv h -v h ®Dr]\dx < e [ rj 2 \Dv h \ 2 dx + [ \v h \ 2 dx . 
Jb+( v ) Jb+( v ) Jb+{ v ) 
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\A(h) ■ Dcp\ dx < e 



Tf \Dv h \ 2 dx+ — 
B+(y) R2 



\v h \ 2 dx 



Bliv) 



(e^tf+L) / (1 + \Du(x + he)\ 2 ) dx , 



B i(v) 



l'B(h) ■ Dip\dx < e 



ce 



r?\Dv h \ 2 dx + ^ 



\vh\ 2 dx 



Bliv) 



+ c([u} c o,x(B+,R")) £ X L 2 



Bliv) 



(1 + \Du(x + he)\ 2 ) dx . 



In order to estimate the last integral on the right-hand side of (|5.2p we calculate 



\ T e-hV\ 

This yields 



'e,-h{ri 2 v h )\ < \h\ 0,3 (\T e ,hu{x - he)\ + \T e , h u(x)\) < 2 [w] C o,x(q+,rw) l/ip 



(5.3) 



/ \h\ 2 \b(x,u,Du) ■ T e ,-h<p\ dx < c([u] C o,a(q+ >r jv)) / (L + L 2 \Du(x)\ 2 ) dx 

JB+(y) JB+iv) 



'Bliv) 

Collecting the estimates for all terms arising in equation 
Caccioppoli-type estimate 



Bt^iv) 



\Dv h ? dx < cFC 2 



B liv) 

L L 2 



\vh\ 2 dx + C 



B liv) 



Bliv) 

and choosing e = ^, we finally conclude the 
(1 + \Du{x)\ 2 + \Du{x + he)\ 2 ) dx , 



and the constant c depends only on ^, and [m]c o a (q+,r n )- With the boundary version of the Sobolev- 
Poincare inequality we deduce 



\Dv h \ dx < c 



B l, 2 iv) 



\Dv h \™+ 2 dx 



Bliv) 



B liv) 



(1 + \Du{x)\ 2 + \Du(x + he)\ 2 ) dx , 



(5.4) 



and the constant c now depends additionally on the dimensions n, N. We here note that the integrand 
of the second integral on the right-hand side of the last inequality belongs to L s °/ 2 due to the higher 
integrability result for Du from (|5.1I) . In the interior we proceed analogously and consider B^(y) with 
centers y G Z {1 _ a)p / 2 (x Q ) satisfying B+(y) C Qt-\h\ and Vn > ^f- If we choose <p := r) 2 (v h - (vh) y ,3R/4,) 
as a test function all the computations above remain valid (with 2 replaced by 4 in inequality (|5.3[) ). Then, 
after applying the Sobolev-Poincare inequality in the interior in the mean value version on the ball B^ji/^{y) 1 
we obtain the corresponding inequality (15.41) with the full ball B R / 2 {y) instead of B^, 2 (y), and c has exactly 
the same dependencies as in the previous reverse Holder-type inequality; in particular, the constant c is 
independent of the parameter h. 

Applying the global Gehring Lemma pTl Theorem 2.4] on the cylinder Zri_ a ) p /2(xo) for the choices of 
a, p and Xq made in the assumptions of the proposition, we find that there exist a constant c depending 
only on n,N,q, ~, [u]c°- x iQ+ .R N ) an d cr and a positive number S depending only on n, N, ~, and 
[u]co,a(q+ rn) such that there holds 



Z„ p ix Q ) 

< c 

< c 



\Dv h \ g dxj " 
(/ \Dv h \ 2 dx 



(l + \Du{x)\ 2 + \Du{x + he)\ 



2\2 



dx 



(l-Ser)p/2 



T e hDu\ dx 



' Z (l~St7)p/lixo) 



(x ) 



Z p /2ixo) 



(1 + \Du(x) 



2\ 2 



dx 
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< c 



(1 + \Du\ 2 ) dxY + [i (1 + \Du\ 2 ) 5 da;) " 



for all q G [2,2 + S). Here, we have also used the bound \h\ < 2ap (with \a\ < jk) and the estimate 
combined with the Holder continuity of u with exponent A (note that as a consequence the constant c then 
depends additionally on the radius p). Hence, for all s £ (2, min{so, 2 + 8}) the previous inequality holds 
true. Keeping in mind the definition of Vh and the higher integrability result (|5.1j) . we finally arrive at 

\T ejh Du\ s dxY <c\h\^(-[ {l + \Du\ 2 )dx 
which finishes the proof of the proposition. □ 



Moreover, we mention two direct consequences of Proposition 15. II The first one follows from Lemma 13.41 
and concerns the slice- wise mean-square deviation of Du: 

Corollary 5.2: Let u £ W^ 2 (Qj,R N ) (1 L co (Q+,R N ) n C Q < X {Q + ,R N ) be a weak solution of flS} under 
the assumptions (H1)-(H3) and (B). Furthermore, let Z p (xq) C Q + for some Xq £ Q + UQ° and a £ (0, j^). 
Then for every 7 £ (0, 1) there exists a function F\ £ L s {Z a p{xo)) such that the following estimate holds 
true: 

\Du(x) — [Du) z i r {x n ) | 2 dx ) 
Z T (z) ' ' ' 

< (/ / \Du{x',x n )-Du(y',x n )\ 2 dy'dx) 2 < cr^ F 1 {z) 

K JZ r {z) JD r (z>) ' 

for all cylinders Z r {z) C Z ap (xo) with z £ Q + U Q°, and the constant c depends only on n, a, A and 7. 

Remark: The L s -norm of F± might blow up if 7 /• 1 (as a consequence of the application of the L q - 
inequality for the maximal operator in the proof of Lemma I3T41 . Moreover, when verifying the assumptions 
of Lemma 13.41 we observe that the number K (resulting from the inequality in Proposition [3TT|) depends on 
the radius p and on a. This dependency is reflected only in the L s -norm of F\. However, this will not be of 
importance because p and a may be chosen fixed in every step of the subsequent iteration. More precisely, 
in the next section we will infer appropriate fractional Sobolev estimates on the cylinders Z ap (xo) and then, 
via a covering argument, also on Q + (respectively on smaller half-cubes in the course of the iteration). 



As a second consequence of Proposition 15.11 we obtain a fractional Sobolev estimate for the tangential 
derivative D'u. This follows immediately from Lemma [3.21 and the inclusion W 0,s C M 8 ' s (for 9 £ (0,1), 
s £ (l,oo)) given in Remarks [ 



Corollary 5.3: Let u £ W^ 2 {QZ,R N ) n L°°(Q£,R N ) n C°> X {Q+,R N ) be a weak solution of (JUT} under 
the assumptions (Hl)-(HS) and (B). Then for every 7 £ (0, 1) there holds 

D'u = (Lqu, . . . , D„_iu) G Af 7QA/2 ^(Q+, M ( "- 1)Ar ) 
for every p < 1. In particular, there exists a function Hi £ L S (Q^ 2 ) such that 

\D'u(x)-D'u{y)\ < \x - y^ (H^x) + H 2 (y)) 

for almost all x, y £ Q^ 2 ' 



5.2 A first estimate for the full derivative 



So far, we can estimate finite differences close to the boundary only with respect to tangential directions. 
In order to find a fractional Sobolev estimate of type (|3.2[) also with respect to the normal direction we next 
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choose a cylinder Z p (xo) C Q + , xq € Q + U Q°, p < p ca cc where p cacc is from Lemma and a £ (0, jq). 
Furthermore, we fix a number 7 € (0, 1) to be specified later. We now study the model system (|4.1I) on 
cylinders Z r {z) with z £ Q + U Q° such that Z% r {z) C Z ap (xo), and by M* we always denote the maximal 
operator restricted to the cylinder Z ap (xo), i.e. 

M*(/)0) := sup / 

for every / £ L 1 (Z ap (xo), R fc ), fc > 1, and z £ Z crp (a;o). We shall frequently use the fact that the maximal 
operator is bounded as a mapping from L p to itself for every p > 1. 



A fractional Sobolev estimate for a n (-,u,Du) 

In coordinates we have the following representation of the weak formulation for the system (|4.1[) : 

N n N 



/Zz~2{ af t (x,u(x),Du(x))D lt <p> dx = ^ f 



b 3 (x, u(x), Du{x)) (p 3 dx 



for all ip £ Cg°(Z r (z), R N ). Following the approach of [35], we are going to derive in the first step a weak 
differentiability result for the function 

A 3 r (x n ) := f a J n (x' ',x n ,u(x' ',x n ),Du(x' ' ,x n ))dx' (5.5) 

JD r {z>) 

for every j £ {1, ...,N} and x n £ I r (z n ). For this purpose we choose a "splitting" test function of the 
form ip(x) — 4>i(x') faixn) Ej where (pi £ Cfi° (D r (z')) with <fii = 1 on the (n — l)-dimensional ball D Tr (z') 
for some r £ (0,1), <f>2 £ (I r (z n )), and where -Ej denotes the standard unit coordinate vector in R^. 
Employing the above identity with such a test function <p then yields 

+ f a j n {x,u{x),Du(x))(j)x{x')D n (f>2{xn)dx' dx n 

Jl r (z„) JD r {z') 

n-1 



1 



^ a j K (x, u(x), Du{x)) D K cj)i{x') <p 2 (xn) dx' dx n 



W I^M*')! J D r (z')\D Tr {z>) K=1 

b 3 (x, u(x), Du(x)) 4>\ (x 1 ) 4>2 (x n ) dx' dx n 



l r (Zn) JD r (z>) 

n-l 



1 



^2 [ a i( x > u(x),Du(x)) - a J K (z, (u) ZiT , (Du) z ,. r (x n ))] 



r(z„) \ D r( z ')\ Jrr JdDf(z') K=1 

x D ti 4> 1 {x')d'K n - 2 {x t )d7(t> 2 {x n )dx n 
b 3 (x, u{x) , Du(x)) 4>\ (x') dx' 4>2 (x„) dx n 



I r (z n ) JD r (z>) 

for j £ {1, . . . , iV}, where we have used the co area formula in the last line. In particular, we may choose 
by approximation a cut-off function of the form 

{1 if \x'-z'\ <rr, 

r-j^f iirr<\x'-z'\<r, 
i{\x'-z'\>r. 

We note that this implies D K cj)i{x') = — ^l T ) r S/Ef?j f° r every n £ {1, . . . , n — 1} provided that \x' — z'\ £ 
(rr,r). Setting 

B 3 K (x) = a j K (x,u(x),Du(x)) - a j K (z, (u) z>r , (Du) z > <r (x n )) (5.6) 
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for j G {1, . . . , N} and k E {1, . . . ,n — 1}, we calculate with this particular choice for cj>\: 

aP n {x, u(x),Du(x)) 4>i (x') dx' D n § 2 (x n ) dx n 
1 



I r (Zn) JD r (z') 



= / \Dl?Y\ f I B ° {x) ■£r J ^dx n - 2 (x')d¥Mxn)dx n 

JI r {z n ) \ U r\Z )\ J tt JdDr(z') \X — Z \ 

+ + f V (x,u(x),Du(x)) (f>i(x') dx' (f>2[x n ) dx n . 

Jl r (z n ) JD r (z') 

Recalling the definition of A 3 r (x n ) given in (|5.5p . we consider the limit t / 1 and conclude from Lebesgue's 
differentiation Theorem that for almost every radius r (and fixed center z G Z <tp (xq)) such that Z r (z) C 
Za-p(xo) there holds 



A 3 r (x n ) D n (j> 2 (x n ) dx r , 



U{z n ) 



\D r (z')\ 



dD r {z') \X — Z \ 



V{X, u{x),Du[x)) dx' 4>2{x n ) dx n ■ 

'l r (z n ) JD T {z>) 

Hence, for almost every radius r with Z r (z) C Z ap (xo) we find that A r (x n ) — (4j(i n ), . . . , A^f (x n )) is 
weakly differentiable on I r {z n ) (note that the index j G {1, . . . , N} and the test function </> 2 are arbitrary 
in the latter identity), and its weak derivative is given by 



A r {Xn) 



1 



\D r (z')\ 



B ( x )--n — ^j7d3{ n - 2 (x') - I b(x,u(x),Du(x))dx' . (5.7) 

> r (z') F _ Z JDr(z') 



ldD r (z') F — * I JD r (z') 

We next consider for any fixed r all radii p G (0, r] and we define the set J via 

2 



J = j p : p G (0, r] and / / 

L Jlo(Zn) Jdl 



Z r {z) 



\B(x)\dx\ 



\B(x)\d'K n -\x')dx n > 

;(z„) JdD^z') r ■ 

The following computations reveal that there holds -C 1 (J) < §: employing the co area formula and Fubini's 
Theorem we get 



/ \B(x)\dx> f f 

JZ r (z) JO Jlp{ 



\B(x)\d3< n -' s (x')dx n dp 



(z n ) JdD ? (z>) 



> 



J Jls(z n ) JdDpiz') 



\B{x)\d'K n ' z {x')dx n dp 



> /-/ \B{x)\dxdp = £}{J)- I \B(x)\dx. 
J.j r J Zr (z) r JzAz) 



Z r {z) 



Therefore, we find some radius p G [5,r] such that on the one hand Ap(x n ) is weakly differentiable and on 
the other hand p £ J ■ Hence, in view of Poincare's inequality and identity (15.71) . we obtain for this choice 
of p: 



{ \Ap(x n ) - (Ap) ZniP \dx n < c(N) \A' p {x n )\dx T 

JlnUru) Jls(gn) 
C(N) 



Ip{z n ) 

< 



\Dp(z')\ Jl p (x n ) JdDp{z') 



5 (z„) 

\B{x)\d'K n ^ 2 {x')dx ri 



c(N) 



Ip(z n ) JDp(z') 



\b(x, u(x), Du{x))\ dx' dx n 



< c{N) 



\Dp{z')\ 



Z r (z) 



B(x)\dx + p+ \b(x, u(x), Du(x))\ dx 

JZp(z) 



< c(n,N) 



f- \B(x)\ dx + r + \b(x,u(x) 7 Du(x))\dx 

Z r (z) JZ r (z) 



(5.8) 
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In the next step we control the integrals arising on the right-hand side of the last inequality by using the 
growth conditions on coefficients and inhomogeneity, respectively, and by exploiting the assumption that u 
is Holder continuous with exponent A (which shall be used without any further comment). 

For the first integral in (|5.8|) we use the definition of B(x) in (|5.6|) . the assumptions (HI), (H3), and 
Corollary [52] to see 



4- \B{x)\ dx < + [ \a(x, u(x), Du(x)) — a(z, (u)z, r , Du{x))\ 

JZ r (z) JZ r (z) 

+ \a(z, (u) z , r ,Du(x)) - a(z, (u) z<r , (Du) z >, r (x n ))\ ] dx 
< AL(r a + Mgo,, (Q+iRN) r aX ) I (l + |I>«|) dx 

JZ r (z) 



L 



Z r (z) 



IZ r (z) 

\Du(x) — {Du) z i _ r {x n ) \ dx 



< cr 2 ^ (M*(l+\Du\)(z)+Fx(z)), 

and the constant c depends only on n, L, [u] c o,\(q+ ,k n ): a i ^ an d 7- Moreover, the functions F\ and M* (l + 
\Du\~) belong to the space L s {Z crp {xo))^ due to Corollary 15.21 and the higher integrability of Du (combined 
with standard properties of the maximal function). 

For the second integral in (|5.8p . we initially assume that we are close to the boundary, meaning that 
z n < 2r. We then infer the following estimate from the natural growth condition (B), the Caccioppoli 
inequality from Lemma 14. II (note that 2r < p ca cc), and the Poincare inequality in the boundary version: 



r+ \b(x,u(x), Du(x))\ dx < r + (L + L2\Du\ 2 ) dx 

JZ r (z) Jz r (z) 



< c[r 



Z r (z) 
1-1 + A 



Du I dx 



Z 2r (z) 



rL 



< cr x M*{\ + \Du\)(z) , 



(5.9) 



and the constant c depends only on n,N, L,L2,v and [u]c o a (q+,k jv )- For cylinders in the interior, where 
z n > 2r, we end up with exactly the same estimate using interior versions of Caccioppoli and Poincare 
where \u\ is replaced by \u — (u) Zt 2r\- 

Hence, combining the last two estimates, we conclude from (|5. 



/ 1/ 

JZp(z n ) Jl 



Dp(z') 



a n (y', x n , u(y', x n ),Du(y', x n )) dy' - 



a n (y,u{y),Du(y))dy 



dx 



= / \A- P {x n ) - (Ap) Zn: p\ dx n < cr 2 ^ [M*(l + \Du\) (z) + F^z)] , (5.10) 

Jl- p (z n ) 

and the constant c depends only on n,N,L,L2,v : [u]c°- x (q+,r n )i a > ^ an d 7- Besides, we have Fx, M*(l + 
\Du\) e L s (Z ap (xo)) for some s > 2. Furthermore, applying Jensen's inequality, conditions (HI), (H3), and 
Corollarv l5.2l we find 



/ 

JZp(z) 



a n (x,u(x), Du(x)) — 



Dp{z') 



a n (y\ x n ,u(y', x n ), Du(y', x n )) dy' 



dx 



< c(L, Mc°.a(q+,r")) p aX I (1 + \F>u\) dx 

JZ- n {z\ 

+ L 1 



Zp(z) 

Du(x' , x n ) — Du{y' , x n ) I dy' dx 

lz p [z) Jd- p (z>) ■ 

< c(n,L, Wco.^Q+^N^a^X^) p 2 ^ [M*(l + \Du\)(z) + F x [z)] 



(5-11) 
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Combining (|5.10|) and (|5.11|) . we conclude 

I \a n (x,u(x),Du(x)) - (a n (- ,u,Du)) \dx < cr 2 ^ \M*(l + \Du\)(z) + Fi(z)] 

for every r with Z r (z) C Z ap {xo) and an appropriate radius p G [§, r] for which Ap{x n ) is weakly differen- 
tiable on I r (z n ) and p J. The constant c here depends only on n, N, L, L 2 , v, [m]c°- a (q+,r jv )j a > ^ an d 7- 
In particular, this yields 

I \a n (x,u(x),Du(x)) - (a n {- ,u,Du)) \dx < cr 1 ^ [M*(l + \Du\)(z) + Fi(z)] , 

JZ r/2 {z) 

and the constant c admits the same dependencies as in the previous inequality. This allows to apply the 
characterization of fractional Sobolev spaces given in Lemma 13.51 and Remarks 13.61 (note that these results 
also hold true if we replace the balls by cubes or cylinders). Since the cylinders Z p {xq) C Q + were chosen 
arbitrarily we infer via a covering argument 

a n (-,u,Du) G M^> S (Q+ 2 ,R N ). 

Furthermore, there exists a function G\ G L s (Q^, 2 ,M. N ) which satisfies 

\a n (x,u(x),Du(x)) - a n (y,u(y),Du(y))\ < \x - yf^" (Gi(sc) + G x {y)) 

for almost every x, y G Qu 2 - We finally note that G\ can be calculated from c, M*(l + \Du\), Fi(z) and 
the restriction on the radius p. 

We close this first step with some remarks concerning the components ak(-,u, Du) of the coefficients, 
fc G {1, . . . , n — 1}, and the interior situation: 

Remarks 5.4: We first note that testing the system (|4.1|) with finite differences in normal direction of the 
weak solution u is not allowed. Hence, the statement in Proposition 15.11 cannot be expected to cover (via 
a modified proof) also differences of Du in any arbitrary direction e G up to the boundary. This 

reveals the crucial point for the up-to-the-boundary estimates derived in this section: the method makes 
only an up to the boundary estimate for a n (-,u, Du) available - which is still sufficient to enable us later 
to find an appropriate fractional Sobolev estimate for Du — but a corresponding estimate for u, Du), 
k G {1, . . . , n — 1}, does not follow. 

For cylinders in the interior, however, Proposition 15.11 holds true for every direction e G As a 

consequence, we may repeat the arguments above line-by-line and end up with an interior fractional estimate 
for the full coefficients a(-,u, Du). We here mention that fractional Sobolev estimates for the coefficients 
a(-,u, Du) are not necessary in the interior to prove the dimension reduction for the singular set. In fact, 
interior fractional Sobolev estimates for weak solutions to elliptic systems with inhomogeneities obeying a 
natural growth condition can be obtained directly by exploiting the fundamental estimate (|4.2j) . see |36j . 

A fractional Sobolev estimate for Du 

The ellipticity condition (H2) and the upper bound in (HI) allow to estimate 

[a n (x,u(x),Du(x)) - a n (x,u(x), Du(y))] ■ (D n u(x) - D n u(y)) 

= J D z a n {x 1 u{x) 1 Du{y)+t{Du{x) ~ Du{y)))dt 

' " (Du(x) - Du{y)) ■ (D n u{x) - D n u(y)) 

> v \D n u(x) - D n u(y)\ 2 - L \D'u(x) - D'u(y)\ \D n u{x) - D n u(y)\ 

for almost all x, y G Q^ 2 ' Dividing by \D n u(x) — D n u(y)\ (provided that D n u(x) ^ D n u(y) which is 
the nontrivial case) and taking into account the fractional Sobolev estimates for both a n (-,u, Du) and the 
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tangential derivative D'u from Corollary 15.31 and condition (H3), the latter inequality implies 

v \D n u(x) - D n u(y)\ < \a n (x,u(x),Du(x)) - a n (x, u(x), Du(y))\ + L \D'u{x) - D'u(y)\ 
<L{\x-y\ a + Mgo,, (Q+)R *) \x - y\ aX ) (l + \Du(y)\) 

+ \x-y\ 2 ^ (Gxix) + Gi(y)) +L\x~y\ 2 ^ (H^x) + H^y)) 
< c(L, Mco,a (q+:R »)) \x - y\^ (1 + \Du(y)\ + G x (x) + G x (y) + H x {x) + H^y)) 



for almost every x, y £ Qt/oi meaning that we have D n u £ Af — - s (Q+ 2 , R N ). Combined with CorollaryO 
we hence end up with 

Du £ M^^{Q+ 2 ,R nN ), 

which is the desired estimate for the full derivative Du. We recall the embedding for the fractional Sobolev 
spaces, namely that 

M^ aX / 2 > s (Q+ 2 ,M nN ) cW^' jaX/2 > s (Q+ /2 ,M nN ) 

for all 7' G (0, 1). Then, in view of the interpolation Theorem 13.31 and the fact that 7 and 7' may be chosen 
arbitrarily close to 1 (an appropriate choice is for example 7 = 7' = ( n +2\ ) 1 ^ 2 )' we finally arrive at the 
higher integrability result 

Du £ L sil+aX/2) (Q+ /2 ,K. nN ). 



5.3 Iteration 

In the next step we iterate the fractional Sobolev estimate for Du. To this aim we define a sequence (bk)keN 
as follows: AAA 

b := 0, b k+1 := ^- + b k (l - - ) = b k + - (a - b k ) 

for all k £ No- We observe that the sequence (b k ) is increasing with b k /• a. The strategy of the proof is 
the following: For every k 6 No we show by induction the following inclusions: 

Du £ L Sk{1+bk \Q+_ k ,R nN ) -> Dti6M*+ 1 ' s '+ 1 (Q+ +i,M nJV ) 

Due L Sh +^ 1+b "+^(Q+ 2k+x ,R nN ), 

where 7 6 (0, 1) is an arbitrary number and where {s k ) k< zfi is a decreasing sequence of higher integrability 
exponents with s k > 2 for every i: £ No. The first step of the induction, k = 0, was already performed 
above (with s± — s). We now proceed to the inductive step: The objective is to find the first inclusion by 
improving the fractional Sobolev estimates in Sections 15.11 and 15.21 and then to deduce in the second step 
the higher integrability result by applying the interpolation Theorem 13.31 

Higher integrability II 

We again need to deduce a higher integrability result for the tangential differences T ei h.Du (cf. Proposi- 
tion l5.1|) which now incorporates the fact that Du is assumed to be integrable with exponent s k (l + b k ). In 
what follows we will frequently use a simple consequence of b k < a, namely the inequality 

a\ + Ml - A) > ^ + b k (l - -) = b k+l . 

Proposition 5.5: Let u £ 2 {Qt , R N ) n L°°{Q+, R N ) n C X (Q+, R N ) be a weak solution to JUJ under 
the assumptions (Hl)-(HS) and (B). Assume further u £ W^ Sk{1+bk) (Q^_ k ,R N ) for some k £ N, s k > 2, 
and let Z p (xq) C Q^-k f or some xq £ Q 2 - k U Q 2 - k , cr £ (0, |), e £ S*™ -1 with e _L e„ and h £ R satisfying 
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\h\ s (0,2(7/o). Then there exists a higher integrability exponent Sk+i G (2, Sk) depending only on n,N, ~, -^7 
and [m]c°. a (<3+,r jv ) such that 



L 



s k+i 

\T eth Du\ Sk+1 dx < c \h\ Sk + lbk + 1 (4 (1 + \Du{x)\Y k(1+bk) dx) 



lZ tTp (xo) JZp(xo) 

for a constant c = c(n,N, £, ^f, [w]c"^(Q+,K^)> P> °) ■ 

Proof: We give only a sketch of proof and refer to [3 Proof of Proposition 8.8] for more details. We start 
from the preliminary estimate (|4.2j) and show that for every 8 G (0, 1) and every cylinder Z r (xo) C Q^-k 
there holds 

f \ Te . h Du\ 2 dx < c\h\ 2bk ^ f (1 + \Du\f +2bh dx (5.12) 

JZ er {xo) Jz r (x„) 

for all e G with e JL e n , /i G R satisfying |/i| < '"C 1 ^" 8 ) and a constant c depending only on 

n, 2V, £, -j^, [u]po,A(Q+ r n), and r. For this purpose, a suitable cut-off function is chosen, and the different 
terms arising on the right-hand side of (I4.2[) are then estimated taking advantage of standard properties of 
finite differences, the integrability of Du with exponent 2 + 2b k and the Holder continuity with exponent A, 
see also [35J p. 387]. 

In the next step we proceed similarly to the case k = and estimate the L Sfc+1 -norm of \T e ^Du\ for some 
exponent Sk+i > 2 in terms of an appropriate power of \h\. To this end we consider directions e G S"™ -1 
with e _L e n and h G R satisfying \h\ < 2~ k ; furthermore, analogously to the proof of Proposition 15 . II we set 

Te,hU ~ -A{h) ~ -3(h) 

and G k (h) = G(h) = fiD z a{x,u(x),Du(x) + tT et hDu(x)) \ dt as above. Analogously to the derivation of 
(15.2(1 we then see that the map Vh,k G W 1 ' 2+2bk (Q^_ k _^,R N ) is a weak solution to a linear system, for 
which the various terms need to be estimated in terms of the I/ 2 -norms of v^.k and Dvh^k- The only point 
differing from the estimates before is the one involving Ti^ih): to find an adequate inequality we first take 
advantage of the Holder continuity of u and Young's inequality and we see 

f (1 + \Du(x + he)\f |r e h u\ 2a dx 
Jb + r ( v ) 

< c([u] c o, HB + r«)) \h\ 2aX ~ 2bkX f (1 + \Du(x + he)\f \T eJl u\ 2bk dx 
Here we have used the fact that 



< c([u] c o, HB+ \h\ 2bk ^ / (I + \Du(x + he)\ + \G h {x)\y + '° k dx . 

JB+{y) 



\r e ,hu\ < \h\ [ \Du(x + the)\ dt =: \h\ G h (x) , 
Jo 

and the function Gh is L Sfc ( 1+b,c )-integrable on B^(y) in view of Fubini's Theorem: 

\G h \ Skil+bk) dx < [ \Du\ Sk( - 1+bk] dx < 00. 
b%(v) JQt- k 

Hence, we find with Young's inequality for every e G (0, 1) 

% {k) {h) ■ DiA dx < L \h\~ bh+1 (l + \Du(x + he)\) \r eh u\ a \Dip\dx 

BJ(b) ' J B+(y) 



<ef vl\Dvi k) \ 2 dx + ^ [ 
JB+ly) X Jb 



lB+(y) JB+(y) 



C 



([u} c o,x {B+RN) )e 1 L 2 / (l + \Du{x + he)\ + \G h (x)\) 2+2bk dx. 
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Arguing exactly as in the proof of Proposition I5.1[ we obtain via the Sobolev-Poincare inequality a reverse 
Holder- type inequality, from which (taking advantage of the higher integrability of Gh) we then deduce the 
desired assertion by the global Gehring Lemma. □ 

Remark 5.6: If Sfe(l + bk) > 2 + 2a is satisfied, it is easy to check that the inequality f|5 . 12[) and in turn 
the statement of the proposition on the higher integrability of the differences hold true with b k +i replaced 
by a. 

Proposition 15.51 combined with Lemma 13.41 and with Lemma 13.21 respectively, again allows to state 
two direct consequences concerning the slice-wise mean-square deviation of Du and a suitable fractional 
differentiability of the tangential derivative D'u: 

Corollary 5.7: Let u G W^(Qt,R N ) D R N ) n C°' A (<2 + , R N ) be a weak solution to $4~B) under 

the assumptions (Hl)-(HS) and (B). Assume further u G W^ Sk{1+bk} (Q+_ k ,R N ) for some k G N, s k > 2, 
and let Z p (xo) C Q 2 -k f or some xq G Q 2 -k U Q^-k and a G (0, g). Then for every 7 G (0, 1) there exists 
a function F k +\ G L Sk+1 (Z ap (xa)) where Sk+i G (2, Sfe) is the higher integrability exponent determined in 
Proposition 15.51 such that the following estimate holds true: 

\Du{x) — (Du) z > r (x„)| 2 da;) 

<(i / \Du(x',x n )-Du(y',x n )\ 2 dy'dxY < cr^"^ F k+1 (z) 

K JZ r {z) JD r {z') J 

for all cylinders Z r (z) C Z ap {xo) with z G Q + U Q°, and the constant c depends only on n, a, A and 7. 

Corollary 5.8: Let u G W^' 2 (Q^,R N ) n L°°(Q%,R N ) n C°> X (Q+,R N ) 6e a weak solution to under 
i/ie assumptions (Hl)-(HS) and (B). Assume further u G Wp^^tQ+^R^) /or some k G N, s fc > 2. 
T/ien /or every 7 G (0, 1) i/iere holds 

D'u G M 7f,fc + 1 ' Sfc + 1 (Q+,R( n - 1 ) Ar ) 
/or every p < 2 ^ +1 . In particular, there exists a function H k +i G L Sk+1 (Qy 2k+1 ) suca that 

\D'u(x) - D'u(y)\ <\x- tf p**+* (H k+1 (x) + H k+1 (y)) 
for almost all x, y G Q^ 2k +i ■ 



An improved fractional Sobolev estimate for a n (-,u,Du) 

Taking into account that Du is assumed to be integrable with exponent s k (l + b k ), we next proceed similarly 
to the case k = 0: We choose a cylinder Z p (xq) C Q^-k with center xo G Qt-k^Q^-k an ^ radius p sufficiently 
small , i. e. p < p ca cc where p caC c is from the Caccioppoli-type inequality in Lemma 14. 1[ and a G (0, |). 
Furthermore, we fix a number 7 G (0, 1) and again study the model system (|4.ip on cylinders Z r (z) with 
z G Q^-fc U Q2- fc sucn that Z2r{z) C Z (Tp {xo). Using the notation from Section [5"T2l we first improve the 
estimate (|5.10[) : To this aim we start with inequality (|5.8[) : For the first integral on the right-hand side of 
(|5.8|) we recall the definition of B(x) in (|5.6|) and take advantage of conditions (HI) and (H3) to infer 

/ \B(x)\dx<L-f (\x - z\ a + \u(x) - (u) z , r \ a ) (1 + \Du{x)\) dx 

JZ r (z) JZ r (z) 




Du(x) — (Du) z i yr (x n ) \ dx . 
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In view of Holder's and Jensen's inequality, the Holder continuity of u and Poincare's Lemma, we derive 
\u{x) - (u) z , r | Q (1 + \Du(x)\) dx 

<(-[ \u(x) - (u) z , r \ a ^ dx)*** ( / (l + |£)u|) 1+i *da:)'* ?FI 

Wz r (z) 1 Wz,.(z) 1 

^ bkX ( I \u(x) - (u) z<r \ 1+b * dx^ll (1 + \Du\) 1+bk dx) **** 



< cr° 



< cf ,aX+b k (l-X) / ^ + 1^1)1+6* dx < cr 7b fc + i M *((! + \Du\) 1+bk )(z) (5.13) 
JZ r (z) 

for c = c(n, [m]c7°- a (q+,r jv ))- Furthermore, we trivially have 

\x - z\ a (1 + \Du{x)\) dx < c{n)r lbk+1 M*((l + \Du\) 1+bk ){z) . 



L 



IZ r (z) 

Keeping in mind Corollary 15. 71 we finally arrive at the following estimate for the integral of |fl(ar)| 



\B(x)\dx < cr 1 '^ 1 (M*((l + \Du\) 1+bk ){z) + F k+1 {z)) (5.14) 

Z,,(z) 

for a constant c depending only on n, L, [v\c°>*(Q+ ,M. N )i ck, A and 7. We mention that the functions M* ((1 + 
\Du\) 1+bk ) and F^+i belong to the space L Sk+1 (Z rTp (xo)), due to the higher integrability of Du and Corol- 
larv 15.71 respectively (note s^+i € (2,Sfc)). 

For the second integral on the right-hand side of (|5.8I) we argue similarly to above on p. 1161 we first 
assume that we are close to the boundary, i. e. z n < 2r. Then, we infer the following estimate from the 
growth condition (B) on the inhomogeneity, the Caccioppoli inequality (note that 2r < p < p ca cc), the 
Holder continuity of u and Poincare's inequality in the boundary version: 



r f \b(x,u(x),Du(x))\dx < r f (L + L 2 \Du\ 2 ) dx 

JZ r {z) JZ r (z) 

<cr ( (l+ ^ 1+ V-M(A-D) d , 
Jz 2r (z) v r ) 

< Cf l+(1-6»)(A-1) f (l+\Du\) 1+bk dx 

JZ 2r (z) 

< cr bk+1 M*((l + \Du\) 1+bk }{z) , (5.15) 

where in the last line we have employed the fact that 1 + (1 — &/.) (A — 1) > bk+i and where the constant c 
depends only on n,N,L,L2,v and [u\c°^(Q+ ,M, N )- For cylinders in the interior, meaning that z n > 2r, we 
end up with exactly the same estimate using both the Caccioppoli inequality and the Poincare inequality 
with \u\ replaced by \u — (u) z ,2r|- 

Merging the estimates found in (15. 14)) and (|5 . 1 5[) together with (|5.8[) hence yields 

JZp(z) JD p (z') 

\A p {x n ) - {A p ) Zn>p \ dx n < cr? 6 ^ 1 [M*((l + \Du\) 1+bk )(z)+F k+1 (z)] 

Ip(z n ) 

for a constant c depending only on n, N, L, L2, [u\c°- x {Q+ ,r n )j a > ^ an d 7- This is the desired improvement 
of inequality (|5TTD]) . Moreover, Ffe+i,M*((l + |Dtt|) 1+6fc ) e L^ 1 (^(xq)) holds true. In order to find 
a fractional Sobolev estimate for the map x 1— > a n (x , u(x) , Du(x)) it still remains to deduce an estimate 



dx 
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corresponding to (|5.11|) . To this aim we follow the line of arguments leading to (|5.11|) and (|5.13j) : in view 
of Corollary 15.71 Holder's inequality and the Holder continuity of u, we see 



f a n (x,u(x),Du(x)) - f a n {y',x n ,u(y',x n ),Du(y',x n ))dy' 



dx 



>D- P {z>) 

<L+ + \Du(x' , x n ) — Du(y' , x n )\ dy' dx 

Jz p (z) Jd- p {z') 

+ 4l/ / (p a + \u(x',x n )-u{y',x n )\ a ) (l+\Du(x)\)dy'dx 

JZp-iz) JD p {z>) 

< cr 7bfc+1 F k+1 {z) + 4Lp a I (1 + \Du(x)\) dx 

Jz- P (z) 

+ 8L(-f \u(x) - (u) z , r \ a1 ^ dx)^ ( / (1 + \Du\) 1+b " dx) 

\Jz r (z) 1 X ~JZ r (z) ' 

< cr 7b, + i (M*((l + \Du\) 1+bk )(z) + F k+1 (z)), 

and the constant c depends only on n, L, [w]c°- a (q+,r jv )i a : ^ ancl 7- I n particular, taking into account 
p G [§? r L we infer from the latter two estimates the inequality 



f \a n (x,u(x),Du(x)) - (a n (- ,u,Du)) ,Jdx 

JZ r/2 (z) 

< cr^ 1 (M*((l + \Du\) 1+b «)(z) + F k+1 (z)) 



where c admits the same dependencies as in the preceding inequalities. Since i^fc+i, M*((l + \Du\) 1+bk ) 
belong to L Sk+1 (Z crp (xo)), we may apply the characterization of fractional Sobolev spaces in Lemma [3751 and 
Remark I3.6[ an d we obtain 

a n ( ■ , u, Du) G M~ /bh+1 ' Sk+1 (Qt /i2 . 2ky R N ) ■ 

Furthermore, there exists G k +i G L Sk+1 (Q^, 2k+1 y R N ) which satisfies for almost every x,y G Q 1 /( 2 fc + 1 ) 

\a n (x,u(x),Du(x)) - a n (y,u{y),Du(y))\ < \x - yp b ^ (G k+1 {x) + G k+1 {y)) . 

We note that G k +\ can be calculated from the constant c, the functions M*((l + \Du\) 1+bk ), ^+1(2) and 
the restriction on the radius p which in turn result in a dependence on the iteration step k. For the interior 
situation we observe that the statements of the Remarks 15.41 remain valid, in particular, the coefficients 
a(-,u,Du) satisfy a corresponding interior fractional Sobolev estimate. 

Final conclusion for Du 

Exactly as before on p. [18] the normal derivative D n u inherits the fractional Sobolev estimate of both the 
coefficients a n (-,u,Du) and the tangential derivative D'u (see Corollary 15. 8|) . This gives 

Du G M^+^+^Q+^^R^). 
At this point we are in the position to use the embeddding 

Ar bk+1 ' Sk+1 {Q+ 2k+1 ,R nN ) c W~<'~ fbk + 1 > Sk + 1 {Q+ 



T®nN\ 
2 k+l j ^ i 



for all 7' G (0, 1). Since 7 and 7' may be chosen arbitrarily close to 1 (the choice 7 = 7' = (^f^x) 1 ^ 2 
like in the first step is still appropriate for every k G N), the application of Theorem 13.31 yields Du G 
LSk+1 (i+b k+1 )(Q+^ k+ ^ R nNy This finishes the iteration and yields: 
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Lemma 5.9: Let a G (0,1) and let u G W^' 2 {Qi,R N ) n L°°(Q%,R N ) n C°' A (Q+, M^), A G (0,1], be a 
weak solution of the Dirichlet problem (|4.1|) under the assumptions (Hl)-(HS) and (B). Then for every 
t < a there exists p — p(t) such that Du G W t ' s (Q^ ,M. nN ), where s > 2 is a higher integrability exponent 
depending only n, N, L, Li, v and [w]c°> x (Q+ ,R N ) no ^ ow 

Proof: In view of bk /* a, the iteration scheme immediately implies a fractional differentiability result 
for Du: for every t\ < a there exists k = k{t\) such that Du G W tl,2 {Qf l ^,R n7V ). In particular, we may 
choose ti sufficiently close to a such that 

2 + 2q< 2 ^ 1 + V. 

In view of Theorem 13.31 we hence find Du G L 2+2a , and at this stage we can indeed stop the iteration: as 
already observed in Remark 15.61 in all calculations above the exponent bj, +1 can be replaced by a, leading 
to 

DueM^iQ+^R^) 

for every t < a and a higher integrability exponent s :— , x > 2. With the choice p := 2~ k ~ 1 the proof of 
the lemma is complete. □ 

Remark: A similar statement was derived for weak solutions to superquadratic nonlinear elliptic systems 
with inhomogeneities satisfying a controllable growth condition, see [18\ Lemma 6.1]. We easily observe 
that the method presented in this section does not only apply to inhomogeneities obeying a natural growth 
condition, but also to those obeying a controllable growth condition. As an advantage of the technique 
presented here, we note that in the formulation of the previous Lemma [53] the low dimensional assumption 
p > n — 2 — 6 for some positive number 8 is not necessary, whereas it was required in the proof of [T8I 
Lemma 6.1]. 

Proof (of Theorem II. ID : Most of the arguments required here can be recovered from [T51 proof of 
Theorem 1.1]; for the sake of completeness we sketch briefly the procedure: First, we reduce the general 
Dirichlet problem (11.11) with boundary values Uq to the corresponding boundary value problem with zero 
boundary values, ie., uq = on <9fi (some attention is needed here: the transformed coefficients still satisfy 
assumption (H1)-(H3), but the transformed inhomogeneity then satisfies also a slightly different critical 
growth condition in the sense that 

\b(x,u,z)\ < L + L 2 (M) \z\ p 

for all (x, u, z) 6 SI x M. N x R nN with \u + uq\ < M. Furthermore, the regularity of dfl allows us to flatten 
the boundary locally around every boundary point xo £ dCl to end up with a finite number of problems of 
type (14.11) on cubes. It then suffices to prove that almost every point on Q° is a regular boundary point, 
i.e. that it belongs to the set Reg Du (Q°): since the Hausdorff dimension is invariant under bi-Lipschitz 
transformations, a standard covering argument then yields the corresponding estimate for the singular 
boundary points on <9S1, i.e. for Sing Du (dfl). 

In the model situation, Theorem 1] guarantees that u is Holder continuous on the regular set 
Reg„(<22 U T) of u with any exponent AG (0, 1 - and that dim M (Sing„(<3^ U Q°) < n - 2. In 

particular, the set of singular points is empty if n = 2. We next observe that the statement in Lemma 15.91 
still holds true if we replace the cube by any smaller cube Q%{xq), meaning that we then obtain 
Du G W t,s (Q~g R (xo), M. nN ) for some S(t) > for alH < a and some s > 2 (independently of the choice of t). 
Therefore, choosing an increasing sequence of sets Bk /" Reg u (Q + U Q°) with Bk C Reg M (Q + U Q°) such 
that Bk is relatively open in Q + U Q° for every k G N, Lemma T5.9I allows us to infer that for every t < a 
and every point x e B k n Q° there holds Du G W^iQ j^Xo), R nN ) for some 5(t) > 0. Taking t G (2/s, a) 
and applying the measure density result [HI Proposition 2.1] (tracing back to Giusti [22l Proposition 2.7]) 
we thus find 

dim M (Sing D JQ°) n Qf(x )) < n - at < n - 1 . 
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Hence, we find dimjf (Smg Du (Q a ) D B^) < n — 2t for every k £ N via a covering argument. Keeping in mind 
dimjf (Sing u ((5 + U Q°) < n — 2, we finally conclude the desired estimate dimj<;(Sing Du ((5 )) < n — 1 on the 
Hausdorff dimension of the singular set for the gradient Du on the boundary. This completes the proof of 
the main result. □ 

Remark: We emphasize that the proof also yields a global fractional differentiability result: if u £ uo + 

is a weak solutions to (jl.lj) under the assumptions of the theorem and with 
a £ (0, 1) arbitrary then Du G W*- 2 (fi, R nN ) for all t < a. 

6 Extensions and open questions related to the dimension reduction 

In the last section we discuss briefly some extensions and open questions related to regularity results 
concerning the above regularity theory and to the Hausdorff dimension of the singular set up to the boundary: 

6.1 The nonquadratic case 

In the case p £ (1, oo) it is convenient to work in terms of the ^-function (useful algebraic properties are 
found in ;13] Lemma 2.1] for the subquadratic case, and for the superquadratic case similar inequalities 
hold true). We start by observing that Mingione [36] was the first to succeed in the dimension reduction 
for the singular points for Du in the interior of O, where u £ W 1,P (H,, W N ) n WL N ) is a weak solution 

(under the standard smallness assumption) to the general system (jl.ll) with an inhomogeneity satisfying a 
natural growth condition and being Holder continuous with respect to every argument: 

\b(x, u, z) — b(x, u, z)\ < L u>p 1 (|x — x\ + \u — u\) (1 + \z\ p ) (6-1) 
\b(x,u,z)-b(x,u,z)\ < L\z-zf 2 (l + |z| 2 + |z| 2 ) p_/32 (6.2) 

for all x,x G u,u £ M. N and z,z £ M. nN , top^t) = minjl,^ 1 } with /?i,/?2 G (0, 1) (these condition are 
in fact only needed to obtain a sharper bound in the superquadratic case). The dimension reduction relies 
on the fundamental estimate (|4.2p which iteratively improves the fractional differentiability of V{Du). In a 
slightly generalized version (including also the subquadratic case) this result can be stated as follows: 

Theorem 6.1 (cf. [36], Theorem 2.2): Let n < p + 2 and u £ W^ p {n, R N ) n L°°(n, R N ) be a bounded 
weak solution to (jl.ip under the conditions (H1)-(H3) and (B), and suppose ||it||L°°(n,R N ) ^ M for some 
M > such that 2 max{p — 1, 1}L2M < v . Then we have 

dimjr ( Sing £ltl (0)) < max{ii - n - 2a, n - p] . 

If p > 2 and a > p/{2p — 2), then under the assumptions (|6.1j) - (|6.2j) for (3\ = 2a — 1 and P2 = (2a — l)/a 
there holds 

dim M (Sing £ , u (fi)) < n - 2a . 

In the general homogeneous case a corresponding interior result [37] is valid, and it turned out that 
the estimates can be extended up to the boundary and that in fact 5{ n_1 -almost every point in fl is a 
regular point for Du independently of the value of p as long as a > 5 and the Holder continuity of u is 
known a priori (J{™ _1 -almost everywhere). For this reason one expects that also in the inhomogeneous case 
under natural growth it should be possible to carry the result of Theorem 16.11 from the interior up to the 
boundary. However, apart from some partial result for p close to 2 (under further restrictions on a), it 
seems impossible to obtain the boundary regularity statement for all p £ (l,oo) with the direct approach 
employed in this paper without any further technical tricks; we emphasize that this problem is not caused 
by the inhomogeneity but appears also for homogeneous systems where boundary regularity was already 
proved by the indirect approach. Hence, the direct approach seems to be matched well only in the quadratic 
situation, and it would be interesting to develop an approach which leads to the known dimension reduction 
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for homogeneous systems, but which can also be applied for inhomogeneous systems under natural growth 
assumption and a general p-growth assumption on the coefficients. 

In fact, the only result available so far in this direction is in the two-dimensional case n = 2, where 
by means of Morrey-type estimates and a comparison principle the optimal (^'"-regularity of every weak 
solution is obtained up to the boundary on the regular set of u (the set where u is locally in a neighborhood 
continuous) , and hence outside a set of Hausdorff dimension less than n — p: 

Theorem 6.2: Let p e (l,oo), Q C R 2 be a domain of class C 1 ^ and u e C 1,Q! (fi, R N ). Let u e 
Mo + Wq' p (£1,M n ) D L°°(n,M. N ) be a bounded weak solution to (p~T|) under the conditions (Hl)-(HS) and 
(B), and suppose that \\u\\lx (n^m 1 *) < M for some M > such that 2max{p — 1, \\L2M < v. Then we 
have 

u E C l ' a (U,M N ) 

in the superquadratic case p > 2, whereas in the subquadratic case 1 < p < 2 there hold 
u e (Reg„ (fi) , R N ) and dim M (ft \ Reg n (H) ) < n - p . 

Proof (Sketch): In the first step of the proof one compares the solution u of the original problem to the 
solution of the frozen, homogeneous problem, for which good a priori estimates are known (in particular 
full C 1 -regularity). This allows to show that Du belongs to the Morrey-space L P,2 ~ T locally on the regular 
set Reg u (ft) of u for every r e (0, 2), see p. 2743] if p > 2 and p. 317] if 1 < p < 2 for the up-to-the- 
boundary versions. This decay estimate for the integrals J(l + \Du\) v dx over balls (possibly intersected 
with SI) in terms of the radius is then used to enter again into the comparison argument and to improve 
the decay of the excess of Du on balls B p (xq) C Br(xq) contained in the regular set Reg u (ft) via 

/ , , \V(Du) - (V(Du)) \ 2 dx 

JB p (xo)nn p ' 

< 2 / \V(Dv)-(V(Dv)) \ 2 dx + 2[ \V(Dv)-V(Du)\ 2 dx, 

where v € u + W 1 ' p {B R (x a ) n fijR^) is the weak solution to diva(x , (u) BR ( Xo - )nn , Dv) — 0. Taking into 
account the Morrey regularity and proceeding as in [331 Chapter 9], we then estimate the first integral on 
the right-hand side via the decay properties for the comparison map (see [TT1 3.12] and [7J (4.34)]) combined 
with a good choice of the radius R as a power of p. The second integral is under control via an estimate 
similar to the ones of (34)] or [8l (4.16)], with the difference that the terms involving s and 6, respectively, 
do not appear if also the higher integrability of the comparison solution is kept in mind. As a consequence, 
the reasoning in |33j applies, and a sufficiently small choice of t shows that the previous excess integral is 
bounded by p 2+a for some a > 0. This corresponds to Holder continuity of V(Du) and therefore of Du 
with some small exponent, which is then improved to the optimal one by standard regularity theory. □ 

6.2 Systems with coefficients a(x, z) 

In the introduction we already spent some words on the situation where the coefficients do not explicitly 
depend on the weak solution u, but only on its gradient and the independent variable. In the interior it 
turned out |37[ I36j that the assumption n < p + 2 (or priori Holder continuity of u) is no longer needed 
(because - roughly speaking - possible singularities of u do propagate to the coefficients only in the gradient 
variable). For this reason the dimension reduction follows in only one step, and it is further valid without 
any restrictions on the space dimension n. In the indirect approach of |18) the same reasoning was applied 
in order to extend these interior results up to the boundary. As a consequence, in case of homogeneous 
systems there holds: 

Theorem 6.3 ([H]): Consider p G (l,oo), a > §. Letfl be a domain of class C 1 '" andu Q E C la (Ti, R N ). 
Assume that u G uq + W ' p (Cl, R w ) is a weak solution to div a(x, Du) = in fl under the assumptions (Hl)- 
(H3). Then < K n ~ 1 -almost every boundary point is a regular point for Du. 
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It is not clear whether the result of Theorem 11.11 can be improved to such vector fields which do not 
explicitly depend on u, in the sense that the existence of regular boundary points in that case is valid for all 
dimensions n > 2. The first problem arises in the preliminary estimate (|4.2p . where in general no positive 
power of h for the last integral - which came from | J b(x, u, Du)T ej -h{r] 2 T ej hu)dx\ - can be produced. 
Hence, to obtain some fractional differentiability of the system, additional regularity assumptions on the 
inhomogeneity are required (such as conditions (|6.1j) - (|6.2[) of Holder continuity with respect to all variables), 
which then allow to use the formula for partial integration for finite differences also in the integral involving 
the right-hand side. In the interior this gives dimj < ;(Sing £)t[ (ri)) < n — 2a in all dimensions. However, in 
the direct approach presented in this paper, these assumptions do not seem to lead to fractional Sobolev 
estimates for x h4 a n {x, Du), see e.g. the derivation of (|5.9[) where the Holder continuity was the crucial 
ingredient. 

6.3 Optimality of the Hausdorff dimension 

It is not clear to what extent the estimates for the Hausdorff dimension of the singular set may still be 
improved, neither in our main statements nor in the cited results [TTJ [SJ [5J ISZJ [3^1 HE] on the Hausdorff 
dimension of the singular set of u and Du. Up to now, the bound for the singular set of Du depends on the 
parameter a. While one cannot rule out that the dependence on a is only due to technique, it is believed 
that this dependence is a structural feature of the problem concerning the Hausdorff dimension of the 
singular set. However, the literature lacks appropriate counterexamples. As a consequence, the question of 
the existence of regular boundary points for Holder exponents a € (0, h) remains open for general nonlinear 
systems of the type considered above. Recently, it was observed by Kristensen and Mingione [33] that the 
Holder continuity assumption in x can be relaxed to a fractional Sobolev dependence. More precisely, still 
assuming the Holder continuity assumption (H3) with now an arbitrarily small a > (as a consequence, 
we have C^'" regularity on the regular set and the characterization of the singular set remains unchanged) , 
we further assume for some /3 > 0: 

(H4) There exists a function g g L v (£l) with g > and v > (1 + /3)//3 such that 

p-i 

\a(x,u, z) - a(x, u, z)\ < L (l + \z\ 2 ) 2 [\x — xf(g(x) + g(x)) + wp{\u — u\)] 

for almost all x,x 6 fl and all u,u g M. N , z 6 R nN (and w^(s) < min{l,s Q }). The function g again plays 
the role of a fractional derivative, see Section [3l Obviously, condition (H4) is weaker than (H3) in the 
case a = ft and v < oo, and it actually turns out that - independently of the value of a - the result of 
Theorem 1 1.1 1 still holds true, provided that j3 > \. In particular, we still get existence of regular boundary 
points (even though in a regular point the exponent of Holder continuity of Du is only a and in general not 

0): 

Theorem 6.4: Consider n G {2,3,4}, a,/3 e (0, 1] with /3 > 1/2. Let tt C R" be a domain of class C 1 ' * 
and u a e C^ Q (r2, R N ). Assume further that u £ u + Wq' 2 (0,,R N ) n L°°{Q,R N ) is a weak solution of the 
Dirichlet problem (|1.1[) under the assumptions (Hl)-(HJ^) and (B), and suppose that \\u\\ Lac ^^N^ < M for 
some M > such that IL^M < v. Then r K a ~ 1 -almost every boundary point is a regular point for Du. 

Proof: The strategy of the proof of the theorem is the same as the one of Theorem ll.il and we immediately 
get into the study of the transformed system (|4.1|) under the assumption of a priori Holder continuity of 
u with exponent A. We first infer Du € L 2+2a {Q^, W lN ) from Lemma 15.91 from a first iteration using 
only assumptions (H1)-(H3) and (B) (alternatively we can use a simple higher intcgrability result via 
Gehring's Lemma), and now start a new iteration as in Section T5.3l bv taking into consideration the additional 
assumption (H4): for this purpose we define a sequence (/3fc)fceN 

and we observe that it is bounded and strictly increasing with limit f3. Observing that only the fractional 
dependence of the coefficients in the x- variable has changed, we now have to re-estimate the terms involving 
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differences of the coefficients with respect to the x- variable, and then all the statement of Section 15731 remain 
true for bk replaced by /3& (and on smaller half-cubes). In fact, there are only two new terms. Under the 
Holder continuity assumption (H3), these were estimated trivially, but they now need to be investigated more 
carefully: the first occurs in the proof of Proposition [53] (assuming that Du is integrable to a power greater 
than 2 4- 2/J&), when we derive a suitable substitute for the preliminary estimate (|4.2|) to find the inequality 
corresponding to (I5.12|) . Actually, only the integral involving A(h) needs to be adjusted: Keeping in mind 
the integrability assumption g e L"(Q^") C 7 '{Qt) an< ^ the Holder continuity of u we calculate with 

(H3), (H4), Young's and Holder's inequality: 

\A(h)\\D{r, 2 T eJl u)\dx < L\hf^^k 

h 

x f (l+|Du(s + fce)|) (g{x + he)+g{x))^ T ^ (r?\T e>h Du\ + 2>q\Dr)\\T e>h u\) dx 

JQ+ 

<e jf\T e:h Du\ 2 dx + <:(£,£) |/i| 2/3 ~w^ +c(L, \\Dti\\l«>) \hf~*~ ^ + 



0+ 



(1 + \Du(x + he)\) 2+2h dx) 1+1 " k ( j \g 
+nspt(>7) ' Wq+ 

where we also have used standard estimates for finite differences (note that the exponent a was for simplicity 
treated as in the powers of \h\). To conclude (|5.12j) it then suffices to observe that both powers of \h\ are 
at least 2/3%+% and that the other terms are estimates exactly as before (but using (H4) instead of (H3) for 
estimating differences of the coefficients with respect to the it- variable) . 

The second new term arises in the fractional Sobolev estimate for a n (x, u, Du) (and in turn in the same 
way also for D n u): it occurs for the first time in the estimate for B(x) and can be dealt with as follows: 



Z r {z) 



\a(x, u{x), Du{x)) — a(z, (u) z>r , Du{x)) \ dx 



<c(r)/ttw[ I (1 + \Du{x + he)\) 1+0k dx + I (g(z) + g(x)) 

'-Jz r (z) JZ r (z) 



1+0 

20 dx 



The right-hand side is then estimated further via the maximal function. 

With these two adjustments, the proof of the theorem then continues as before, leading to the existence 
of regular boundary points for a £ (0, 1) arbitrarily, provided that /? > 4. □ 
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